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ABSTRACT 



Context. Several recent spectrophotometric studies failed to detect significant global magnetic fields in late-B HgMn chemically pe- 
culiar stars, but some investigations have suggested the presence of strong unstructured or tangled fields in these objects. 
Aims. We used detailed spectrum synthesis analysis to search for evidence of tangled magnetic fields in high-quality observed spectra 
of 8 slowly rotating HgMn stars and one normal late-B star. We also evaluated recent sporadic detections of weak longitudinal mag- 
netic fields in HgMn stars based on the moment technique. 

Methods. Our spectrum synthesis code calculated the Zeeman broadening of metal lines in HARPS spectra, assuming an unstruc- 
tured, turbulent magnetic field. A simple line formation model with a homogeneous radial field distribution was applied to assess 
compatibility between previous longitudinal field measurements and the observed mean circular polarization signatures. 
Results. Our analysis of the Zeeman broadening of magnetically sensitive spectral lines reveals no evidence of tangled magnetic fields 
in any of the studied HgMn or normal stars. We infer upper limits of 200-700 G for the mean magnetic field modulus - much smaller 
than the field strengths implied by studies based on differential magnetic line intensification and quadratic field diagnostics. The new 
HARPSpol longitudinal field measurements for the extreme HgMn star HD 65949 and the normal late-B star 21 Peg are consistent 
with zero at a precision of 3-6 G. Re-analysis of our Stokes V spectra of the spotted HgMn star HD 11753 shows that the recent 
moment technique measurements retrieved from the same data are incompatible with the lack of circular polarization signatures in 
the spectrum of this star 

Conclusions. We conclude that there is no evidence for substantial tangled magnetic fields on the surfaces of studied HgMn stars. We 
cannot independently confirm the presence of very strong quadratic or marginal longitudinal fields for these stars, so results from the 
moment technique are likely to be spurious. 

Key words, stars: atmospheres - stars: chemically pecuhar - stars: individual: HD 35548, HD65949, HD71066, HD78316, 
HD 175640, HD 178065, HD 186122, HD 193452, HD 209459 - stais: magnetic field - polarization 
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Mercury-manganese (HgMn) stars comprise a group of late- 
B chemically peculiar stars distinguished by a strong over- 
abundan ce and an unu s ual isotopic composition of he avy el- 
ements t Adelman etall l2004t IWoolf & Lambertl [1999). Many 
of these stars are slow ro tators and members of close spec- 
trosco pic binary systems (lAbt et al.L l2002t ICatanzaro & LetoL 
12004 !). The properties of HgMn stars facilitate precise chemi- 
cal a bundance analysis of their atmospheres (e.g. Adelman et al., 
12006 ). making them the preferred targets for detailed com- 
parisons of observed surface abundan ce patterns and pre- 
dictions from atoniic diff usion theory (iMichaud et al.L Il974t 
lAlecian & MichaudLll98lh . 

The unexpected discovery of low-contrast abundance inho- 
mogeneities on the surfaces of HgMn stars has rekindled in- 
terest in the ast r ophys i cal processes that operate in these stars 
(Adelm an etall 120021: iKochukhov et all 120051: iHubrig et all 
r2006at IKochukhov et al.L l2011ah . Chemical spots are often 
found in magnetic B-type (Bp) stars, which overlap with HgMn 



* Based on observations collected at the European Southern 
Observatory, Chile (ESO programmes 084.D-0338, 085.D-0296, 
086.D-0240). 



Stars on the H-R diagram. Bp stars possess strong global mag- 
netic fields, which are believed to be responsible for the chem- 
ical spot formation (e.g. Michaudet al., 1981). In contrast to 
this clear link between magnetic fields and chemical spots 
in Bp stars, no robust and reproducible magnetic field de- 
tections have ever been reported for any of the HgMn stars. 
Moreover, inhomogeneities on HgMn stars evolve with time 
(Kochukhov et al., 2007; Briquet et al., 2010), possibly indicat- 
ing a previously unkno wn time-depen dent, non-equilibrium dif- 
fusion process (A lecian et al.L 1201 lb . This behavior is, again, 
very different from that of spots on magnetic Bp stars, which are 
obser ved to be stable for at least several decades (lAdelman et al.L 
|200T1) . 

The possible role of weak magnetic fields in explaining 
the puzzling surface phenomena observed in HgMn stars is 
a topic of ongoing debate. The literature contains claims of 
magnetic field detections based on low- and high-resolution 
spectr opolarimetric observations (iHubrig et al.l l2006bl 1201 Ol 
l2012h . However, many other detailed circular polarization 
studies of individual spotted HgMn stars revealed n o mag- 
netic field si gnatur e s in s pectral line profiles (Wade et al. , 
2006; Folso metaLL l2()10l; iMakaganiuk et all 1201 lal 12012 : 
IKochukhov et all|201 lah . Similar null results were also reported 
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by several spectropolarimet ric surveys that included a large 
number of HgMn-type stars ("Sh orlin et al.l l2002t lAuriere et all 
[2010; Makaganiuket al., 2011b). These studies established an 
upper limit from a few tens of G to just a few G for the mean 
longitudinal magnetic field in HgMn stars. 

The failure of spectropolarimetric studies to unambiguously 
detect the fields in HgMn stars may be attributed to t he com- 
plexit y of the surface magnetic field topologies (e.g. iHubrigL 
Il998h . Indeed, spectropolarimetry and, especially, the mean lon- 
gitudinal magnetic field diagnostic is only sensitive to the line- 
of-sight magnetic field component. Theoretically, the net line- 
of-sight magnetic field can be close to zero if many regions 
of opposite polarity are present on the stellar surface. Highly 
structured magnetic field topologies can reduce the Stokes V 
signatures in spectral line profiles below the detection thresh- 
old. On the other hand, field orientation has a relatively minor 
effect on magnetic broadening and splitting, so this diagnos- 
tic can reveal complex magnetic fields invisible to spectropo- 
larimetry. This type of analysis has yielded a number of surpris- 
ingly high magnetic field estimates for HgMn stars. For instance, 
on the basis of studying the differential magnetic intensification 
of Fe 11 lines, iHubrig et al.l (1 19991) and lHubrig & Castellil (1200 lb 
concluded that some HgMn stars possess complex or "tangled" 
magnetic fields stronger than 2 kG. A complementary quadratic 
magnetic field diagnostic method based on the comparison of 
magnetic broadening in spectral lines with different Zeeman 
splitting patterns suggested the prese nce of even stronger, 2- 
8 kG, magnetic fie lds in HgMn stars (iMathvs & Hubrigj [T995t 
lHubrigetal.ll2012h . 

Both the magnetic intensification and quadratic field diag- 
nostic methods rely on a number of simplifying assumptions 
about the spectral line formation in magnetic field. These meth- 
ods were originally developed for intensity spectra with lim- 
ited wavelength coverage and for spectropolarimetric observa- 
tions with moderate resolution and relatively low signal-to-noise 
(S/N) ratio. Modem echelle spectra of HgMn stars provide ac- 
cess to numerous magnetically sensitive lines that can be ana- 
lyzed using sophisticated theoretical tools. The goal of our paper 
is to probe the existence of complex magnetic fields in HgMn 
stars by taking advantage of the high-resolution, high S/N ob- 
servations and using accurate methods to model magnetic ra- 
diative transfer We also critically examine spectropolarimetric 
magnetic field detections that seem to contradict our null results 
for a few HgMn stars (Makaganiuk et al., 201 lab, 2012). 

The rest of the paper is organized as follows. Sect. |2]presents 
our spectroscopic and spectropolarimetric observations and dis- 
cusses corresponding analysis methods. Results of the magnetic 
field search using circular polarization measurements, differen- 
tial magnetic intensification, and Zeeman broadening of spectral 
lines are reported in Sect.|3] Finally, Sect. |4]compares our results 
with the outcome of previous investigations and assesses impli- 
cations of our study for the general question of the presence of 
magnetic fields in HgMn stars. 



2. Methods 

2.1. Observations and data reduction 

We analyze spect ra of HgMn stars obtained with the circu lar po 



(l2011bllaLl2012h . We also obtained new circular polarization ob- 
servations of the HgMn star HD 65949 and analyzed spectra of 
the normal late-B star HD 209459, which was observed together 
with the HgMn st ars from our sample but was not studied by 
iMakaganiuk et all 

The H ARPSpol instrument records spectra in the wavelength 
range 3780-6913 A with an 80 A gap around 5300 A. All stars 
in our study were observed using the circular polarization an- 
alyzer. Observations of each star were typically split into four 
sub-exposures between which the quarter-wave retarder plate 
was rotated in 90° steps relative to the beamsplitter The spec- 
tra were extracted a nd cahbrated using a dedic ated version of 
the REDUCE pipeline dPiskunov & ValentiL l2002l) . The Stokes / 
spectrum was obtained by averaging the right- and left-hand po- 
larized spectra from all sub-exposures of a given star. The Stokes 
V and diagnostic null spectra were derived by combining the ex- 
tracted spectra accor ding to the "ratio" spectropolarimetric de - 
modulation method (iDonati et all 119971: iBagnulo et all l2009l) . 
Analysis of the emission lines in the ThAr comparison spec- 
tra showed that the resolving power of HARPS in polarimetric 
mode is /? = A/AA - 109 000 with about 1-2% variation across 
the echelle format. 

Subtle Zeeman broadening due to a weak and complex mag- 
netic field can be reliably detected only for very slowly ro- 
tating stars. Therefore, we selected a small number of sharp- 
lined targets out of the full sample of HgM n stars examined 
in the survey by IMakaganiuk et al.l (l2011bl) . Our primary tar- 
gets include 7 HgMn stars with VeSin/<2 kms"': HD 35548, 
HD 65949, HD 71066, HD 175640, HD 178065, HD 186122, 
and HD 193452. For comparison, we also analyzed two stars 
with projected rotational velocities in the range of 4-6 kms"': 
HD78316(;^Cnc) and HD 209459 (21 Peg). Table[I] summarizes 
alternative designations, atmospheric parameters and literature 
information on the rotational velocities of our program stars. 

Preliminary continuum normalization of Stokes / spectra 
was performed by consecutively applying the blaze function, an 
empirical response function, and sm ooth fit corrections as de- 
scribed bv IMakaganiuk et all (1201 Ibh . The resulting continuum 
normalization has a precision of approximately 0.5%. The final 
zeroth-order correction to the continuum was applied as part of 
the modeling of intensity profiles of individual spectral lines de- 
scribed below (Sect. l331 l. 

Characteristics of the HARPSpol observations of 9 program 
stars are provided in Table |2] We list the UT dates of individ- 
ual observations, the corresponding heliocentric Julian dates, 
the total exposure time, and the peak S/N per pixel measured 
at /i x; 5200 A. Two observations were obtained for 21 Peg. 
Both were used to measure the mean longitudinal magnetic field 
(Sect. ITT), but only the higher quality one was employed for the 
Stokes / Zeeman broadening analysis (Sect. l33T l. 



larimetric mode (ISnik et all 20111 



HARPS spectrometer ("Mav or et al. 



Piskunovetail 1201 lb of the 
l2003l) installed at the 3.6-m 



ESO telescope in La Silla. Most observations were obtained dur- 
ing several observing runs in 2010. The acqui sition and reduc- 
tion of these data was described in detail by Makaganiuk et al.l 



2.2. Muiti-line polarization analysis 

A search for weak magnetic field signatures in high-resolution 
stellar Stokes V spectra is greatly facilitated by a multi-line anal- 
ysis. In this study we used the least-squares deconvolution (LSD) 
code developed by Kochukhov et al. (2010). The LSD method, 
originally introduced bv .Donati et al..(.1997.) . combines informa- 
tion from all suitable stellar absorption lines under the assump- 
tion of self-similarity of the line profile shapes and a linear addi- 
tion of overlapping lines. Kochukhov et al. (2010) showed that 
these assumptions are adequate for weak and moderately strong 
magnetic fields and that the mean longitudinal magnetic field 
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inferred from the Stokes V LSD profiles is not affected by sys- 
tematic biases. 

In this study we applied the LSD analysis to our new spec- 
tropolarimetric observations of HD 65949 and HD 209459. For 
the remaining targets the LSD longitudinal field measurements 
were published by Makaganiuket al. (2011b). For complete- 
ness, Table|2]reproduces these measurements. The line lists nec- 
essary for app l ying L SD were obtained from the VALD database 
(JKupka et all Il999h . using atmospheric parameters listed in 
Table [1] and adopting chemical abundances from the studies by 
[Cowley et al. (2010) and Fossati et al. (2009) for HD 65949 and 
HD 209459, respectively. In total, we used 500-600 lines deeper 
than 10% of the continuum for each star. The application of LSD 
increased the S/N of the Stokes V profiles by a factor of ^ 6. 
Sect.[TT]describes the analysis of the LSD profiles of HD 65949 
and HD 209459. 



2.3. Magnetic intensification and broadening of spectral lines 

The Zeeman effect produces intensification, broadening, and 
splitting of a spectral line depending on its effective Lande fac- 
tor and the magnetic splitting pattern. However, unless the field 
is very strong and the Zeeman splitting can be resolved and mea- 
sured directly, there are few, if any, reliable model-independent 
methods to derive magnetic field intensity from the Stokes / 
spectra. In particular, as we will show below, such techniques as 
differential line intensification and quadratic field diagnostic are 
prone to biases and lead to contradictory results for late-B stars. 
In this situation, detailed spectrum synthesis is the only robust 
way to dete ct and me asure mag netic fields with strengths be - 
low 1-2 kG (Kochukhov et al., 2002, 2006; Johns-Krull. '20071) . 
Refined analysis of high-resolution spectra of bright stars can 
even reveal fields in the 400-500 G range (lAnderson et al.l 
I2010I) . 

Here we use the syn mast magnetic spectrum synthesis code 
dKochukhov et 311120101) to simulate the effects of different mag- 
netic field geometries on the profiles of magnetically sensitive 
metal lines. This code calculates the four Stokes parameter spec- 
tra for a homogeneous magnetic field distribution characterized 
by the three vector components specified in the stellar coordinate 
system (see JKochukhovl 1 20071) or for a turbulent magnetic field 
as described below. These calculations are based on the llmod- 
ELS model atmospheres (Shulvak et al., 2004), computed for the 
stellar parameters given in Table [T] taking into account individ- 
ual abundances of the program stars. We set the micro turbulent 
velocity to zero, since HgMn stars are known to lack sign atures 
of turbulence in their atmospheres (JLandstreet et al.l l2009l) . 

As mentioned in the introduction, previous studies of HgMn 
stars could not unambiguously detect Zeeman-induced polar- 
ization in spectral lines. Even if we believe a few disputed de- 
tections, the polarimetrically-inferred magnetic fields in HgMn 
stars typically do not exceed ~ 100 G. At the same time, multi- 
kG fields were inferred from the analysis of intensity spectra of 
the same or similar stars. These vastly different field strength es- 
timates can be reconciled only if the fields are arranged on the 
stellar surfaces in a highly structured, complex configuration. To 
model such fiel ds we use the concept of an isotropic, tu rbulent 
magnetic field (JLandi DeglTnnocenti & Landolfil l2004l) . which 
implies magnetic field parameters that are uncorrelated on the 
characteristic scale of radiative transfer A turbulent field repre- 
sents an opposite extreme compared to an organized, large-scale 
magnetic field topology. Yet, it is the only known magnetic field 
model satisfying the observational requirement of giving a no- 



Fe II 4273.3 A, Vsini= 2 km/s 



Fe II 6149.3 A, Vsini= 2 km/s 




A (A) 



Fe II 4273.3 A, Vsini=10 km/s 



A (A) 
Fe II 6149.3 A, Vsini=10 km/s 






1.4 








1.2 




„„,,,„,,„- 


& 


1.0 




^^ 


B 
c 


""'^^'■^.s^ 




0.8 




- 






B =0.0 kG 


B=1.5kG 




0.6 


B =0.5 kG 
B=1.0kG 


B, = 2.0kG 
B=3.0kG 



4273.0 4273.5 

A (A) 



6149.0 6149.5 

A (A) 



Fig. 1. Effect of the turbulent magnetic field on the Fen 4273.3 
and 6149.3 A lines. The solid curves show synthetic profiles for 
the field strengths ranging from to 3 kG and for Vg sin / = 2 
and 10 km s ' . The dotted lines show the difference between the 
magnetic and non-magnetic profiles, offset vertically by 1 .2. The 
error bar in the upper right corner of each panel corresponds to 
the S/N ratio of 200:1. 



ticeable Zeeman broadening and intensification in the intensity 
spectra without producing a detectable circular polarization. 

In general, the isotropic, turbulent magnetic field is char- 
acterized by some distribution of the magnetic field modulus, 
F{Bi). For this magnetic field topology, the polarized radiative 
transfer equation terms that depend on the field orientation av- 
erage to zero, yielding null emergent QUV spectra but retain- 
ing the magnetic field effects in Stokes /. We implemented this 
model in the synmast code assuming, as a first approximation, 
that the entire stellar surface is covered by the magnetic field 
with a strength {B^}. The disk-integrated spectra for arbitrary 
Ve sin / are then produced from the intensity calculations at sev- 
eral limb angles, similar to a non-magnetic spectrum synthesis 
(Kochukhov. 2007) . 

Figure [Uillustrates the effect of a turbulent magnetic field on 
the Fell 4273.3 (z = 2.15, pseudo-triplet spHtting) and 6149.3 A 
(z = 1.35, doublet splitting) spectral lines. These calculations 
were done for the T^ff - 12000 K, log g - 4.0 model atmosphere 
and a spectral resolution R = 109 000. As one can see from this 
figure, for a very low Ve sin /, both lines show detectable profile 
distortions for (B,)> 0.5-1.0 kG fields. Starting from ~2 kG, 
the partially-resolved Zeeman splitting becomes obvious in the 
Fell 6149.3 A line. For higher projected rotational velocities a 
stronger field is required to produce significant line profile dis- 
tortions. Based on these calculations, we conclude that, provided 
Ve sin / can be independently estimated from magnetically insen- 
sitive lines, it is possible to detect kG-strength tangled magnetic 
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fields using magnetically sensitive spectral lines in S/N > 200, 
R > 10^ observations of slowly rotating late-B stars. 



2.4. Modeling of circular polarization profiles 

Historically, the efforts to detect magnetic fields in early-type 
chemically peculiar stars mostly refied on the l ongitudinal 
field diagnostic (e.g. lBorra & LandstreeA Il980l iMathvs. .1991.) . 
However, modem availability of the high-quality mean Stokes 
V profiles enables several other field detection and assessment 
methods. For instance. Petit & Wade (2012) demonstrated how 
a useful estimate of the surface dipolar field strength can be ob- 
tained from the LSD Stokes V profiles w ithout complete rota- 
tional phase coverage. IShultz et aP (l2012h used a similar polar- 
ization modeling approach to test compatibility of their high- 
resolution Stokes V observations with the magnetic field mod- 
els derived solely from the (B^) measurements. Here we apply 
a similar technique to assess the circular polarization signatures 
expected for individual (B^) detections reported in the literature. 

To model the observed LSD Stokes / and V profiles, we 
use the polar i zation spectrum synthesis method described by 
iPetit & Wadd (l2012h . The local intensity profile is represented 
by a Gaussian function while the Stokes V is given by the scaled 
derivative of the Gaussian ac cording to the well-known weak- 
field approximation formula (ILandi Degl'Innocenti & LandolfiL 
12004) . The wavelength and efi^ective Lande factor required by 
this approximation are given by the respective mean values for 
the LSD mask. The Stokes / and V spectra are obtained by nu- 
merical integration of these local profiles over the visible stellar 
hemisphere for a given magnetic field geometry. For the calcula- 
tions presented in this paper, we adopted a homogeneous radial 
field distribution with the strength adjusted to match a particular 
(Bz) value. Other free parameters of the model (vg sin / and line 
depth) were obtained by fitting the observed Stokes / profile. 

Although not as realistic as e.g. dipolar field, such schematic 
field model has the advantage of producing the same basic anti- 
symmetric Stokes V profile shape independently of the viewing 
angle. Experimenting with different magnetic field distributions, 
we found that a uniform radial field provides the most conser- 
vative estimate of the circular polarization amplitude for a given 
(Bz). In other words, any other surface magnetic field distribu- 
tion must produce a stronger Stokes V signature for the same 
(Bz)- This is illustrated in Fig. |2] where we show the Stokes V 
profiles computed for a perpendicular dipolar field and for radial 
fields with the same (B^) at each rotational phase. Evidently, the 
radial field distributions yield a lower amplitude of the circular 
polarization signal. This difference is small for very low Ve sin / 
but then quickly increases for larger projected rotational veloci- 
ties because of the crossover effect ( Mathvs..l995a) . Thus, one 
can ascertain that the Stokes V profiles corresponding to a homo- 
geneous radial field provide a robust lower limit on the circular 
polarization signal that must accompany a given {B^} measure- 
ment. 



3. Results 

3.1. Non-detection of the longitudinal magnetic field in 
HD 65949 and HD 209459 

The outcome of our single Stokes V observatio n of the ex- 
treme HgMn star HD 65949 (ICowlev et all 12010") is illustrated 
in Fig. |3] The LSD Stokes V profile has no polarization signa- 
ture stronger than 10"^. The mean longitudinal magnetic field, 
{Z?z) - -2.8 + 4.6 G, inferred from the first moment of the cir- 
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Fig. 2. Comparison of the Stokes V signatures for the dipolar 
field topology with i - p - 90° {solid lines) and for a homoge- 
neous radial field distribution (dashed lines) scaled to give the 
same mean longitudinal field at each rotational phase. The left 
and right panels show results for v'eSin/= 1 and 10 kms"', re- 
spectively. The Stokes V spectra are offset vertically according 
to the stellar rotational phase. 



cular polarization profile is also formally consistent with zero. 
Thus, we conclude that there is no evidence of a magnetic field 
on HD 65949. 

iHubrig et al.] (1201 2h published five FORSl/2 longitudinal 
magnetic field measurements of HD 65949 with typical error 
bars of 20-60 G. Four of their (Bz) values, distributed over 7 
years, sugges t a magne t ic field detection in the range -77 G to 
-290 G. ,BagnuloetalJ(l2012h reanalyzed the FORSl observa- 
tion from which the last (and the strongest) of these measure- 
ments was derived, finding (Bz) - -42 + 79 G. They conclude 
that in this and many other cases, the FORS detections of weak 
magnetic fields are not reliable. 

We used the method described in Sect. 12.41 to predict the 
St okes V signatures co rresponding to the (Bz) values reported 
by Hubr ig et al.l (|2012'). As demonstrated by the dashed lines 
in Fig. |3] a longitudinal field at the level of -77 to -290 G 
should produce a prominent polarization signature, readily de- 
tectable with the LSD profile quality achieved in our observa- 
tions. It is highly improbable that this huge difference between 
the predicted and observed polarization signature is due to ro- 
tational modulation. Thus, the HARPSpol data are incompatible 
with previous claims of the longitudinal magnetic field detec- 
tions in HD 65949. 

Two HARPSpol spectropolarimetric observations are avail- 
able for the superficially normal late-B star HD 209459 (21 Peg). 
The LSD Stokes / and V profiles derived from these spectra are 
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Table 1. Atmospheric parameters and rotational velocities of HgMn stars. 



HD 



35548 
65949 
71066 
78316 



186122 
193452 
209459 



Star 
HR 



Other 



T'eff log g Reference v^ sin ; Reference v^ sin i, this study 
(K) (cgs) (kms-') (kms-') 



1800 

3302 
3623 



y^Vol 
kChc 



175640 7143 
178065 7245 



7493 46Aql 
7775 f3~ Cap 
8404 21 Peg 



11100 
13100 
12000 
13200 
12000 
12300 
12560 
10750 
10400 



3.80 
4.00 
4.10 
3.70 
3.95 
3.60 
3.80 
4.00 
3.55 



1 

3 

4 

10 

5 

6 

7 

12 

9 



1.0-2.0 
0.0-1.0 
1.5-2.0 
6.5-7.3 
1.5-2.5 
1.5-2.0 
0.0-2.0 
0.75-2.0 
3.7-3.8 



1,2 

3 

2,4 

1,6,11 

1, 5, 6, 8 

1, 2, 6, 8 

1, 2, 6, 7, 8 

1,2, 

6, 



,6,8 
,8,9 



1.0 ±0.2 
1.3 ±0.3 
1.7 ±0.2 

6.5 ±0.3 

1.6 ±0.3 
1.9 ±0.2 

1.1 ±0.2 

1.2 ±0.2 
3.6 ± 0.2 



References. (1) [Polk et alljlooa), (2) iHubrig et al.l ( fT99^, (3) ICo wlev et al.l . , . . - . - . -- - ,, . , 

Lands treet et al.1 KOOgh. (7) ICastelli et al.' ('2009'). (S) iHubrig &Castelli. (.200 ih . (9) Ipossati et al.l (l2009h . CIO) iRvabchikova et al.. a998.^ . (11) 
Wool f & Lambert! J1999h , (12)|Wahleren et al. (200(|. 



■ ( l2Ol0h, (4)lYuce et al.l dlQUl) . ( 5) ICastelli & Hu brig ("200?), (6) 



Table 2. Observational data and results of the magnetic field analysis of HgMn stars. 



Star 



UT Date 



HJD -24 X 10' 



rexE_(s) 



S/N {B,} (G)" (fit) (G)" 



HD 35548 


2010-01-06 


55202.7401 


920 


200 


-0.7 ± 2.5 


<200 


HD 65949 


2011-02-08 


55600.7080 


2400 


180 


-2.8 ±4.6 


<300 


HD 71066 


2010-01-09 


55205.7867 


2000 


400 


-1.1 ±0.8 


<200 


HD78316 


2010-01-15 


55211.7873 


1200 


370 


-0.4 ± 2.7 


<700 


HD 175640 


2010-05-01 


55317.8656 


1200 


200 


-0.6 ± 2.2 


<250 


HD 178065 


2010-05-03 


55319.8317 


1600 


220 


-2.3 ±2.1 


<500 


HD 186122 


2010-05-03 


55319.8505 


1400 


230 


0.7 ± 1.9 


<200 


HD 193452 


2010-05-03 


55319.8683 


1200 


240 


-1.0± 1.4 


<250 


HD 209459 


2010-08-09 


55417.6295 


1600 


90 


-2.2 ± 6.0 






2010-08-13 


55421.7623 


800 


170 


-4.6 ±3.1 


<600 



Notes. '"' Longitudinal magnetic field from iMakaganiuk et al] ( 1201 Ibl) except for HD 65949 and HD 209459. '" Upper limit on the turbulent 
magnetic field strength at the 95% confidence level. 
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Fig. 3. LSD Stokes / and V profiles for HD 65949. The mean po- 
larization profile is offset vertically and expanded by a factor of 
10 relative to Stokes /. The dashed lines show the circular polar- 
ization signatures for (B^) = -77 and -290 G, coiTesponding to 
th e range of th e FORSl/2 longitudinal field detections reported 
bv lHubrig et~al. t2012i) . 
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Fig. 4. LSD Stokes / and V profiles for the two HARPSpol ob- 
servations of HD 209459. The mean polarization profiles are off- 
set vertically and expanded by a factor of 10 relative to Stokes 
/. The dashed lines show the circular polarization signatures for 
(Bz) = 53 and 99 G, coiTespondi ng to the range of t he longitu- 
dinal field detections claimed by IHubrig et al.l (1201 2i) from the 
same observational data. 



presented in Fig. |4] No evidence of the Zeeman Stokes V sig- 
natures can be found above the noise level of 1.7 x 10"^ for the 
first (9 Aug 2010) and 8.7 x lO"'* for the second (13 Aug 2010) 
observation. The corresponding longitudinal magnetic field is 



-2.2 + 6.0 G and -4.6 + 3.1 G, respectively. Thus, it is unlikely 
that 21 Peg possesses a large-sca le field stronger than a few G. 

Contrary to these results, Hu brig et al.l (1201 2b reported a 50- 
100 G longitudinal field based on the moment technique anal- 
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ysis (iMathvsL [19911) of the same HARPSpol spectra of 21 Peg. 
Magnetic field was presumably found in the lines of Ti, Cr, and 
Fe, which comprise 80% of the features included in our LSD 
line mask for this star. Surprisingly, these detections were re- 
ported only for the lower quality spectrum but not for the higher 
quality one obtained four days later The minimum Stokes V 
signature corresponding to the {B^) measurements reported by 
[Hubrig et al. (2012) significantly exceeds the noise level of the 
LSD Stokes V profiles (see Fig. |4]i. The lower S/N observa- 
tion is incompatible with the predicted polarization signature for 
(Bz) = 53 G at the confidence level of » 99.9% according to a.x^ 
probability analysis. The lack of polarization signal in the data is 
therefore grossly inconsistent with the longitudinal field detec- 
tions obtained for 21 Peg using the moment technique. 

3.2. Magnetic intensification of the Fe ii 6147 and 6149 A 
spectral lines 

The singly ionized iron absorption lines at 6147.7 and 6149.3 A 
share the same upper atomic energy level and have nearly identi- 
cal oscillator strengths. However, their Zeeman splitting patterns 
are very diff'erent: the 6147.7 A line splits in two n and four cr 
components whereas the 6149.3 A line has two cr components 
coinciding in wavelength with the pair of tt components. As a 
result, the latter line exhibits a simple doublet Zeeman splitting 
pattern, which makes this line particularly useful for direct mea- 
surement s of the mean m agnetic field modulus in slowly rotating 
Ap stars (Ma thvs et al.L Il997h . despite its relatively small effec- 
tive Lande factor (z = 1 .35). 

The two Fe ii lines have very similar intensity in normal stars. 
On the other hand, in the presence of the magnetic field with a 
strength above a few kG, the Feii 6147.7 A line experiences a 
stronger magn etic i ntensification compare d to the 6149.3 A line. 
iMathvsl (Il990h and lMathvs & Lanzl ( 1 19921) noted that in the field 
strength interval from 3 to 5 kG the normalized equivalent width 
difference of these lines, AWa/{Wa} = 2(W6i47-W6i49)/(W6i47 + 
^6149), follows an approximately hnear tre nd with mean mag- 
netic field modulus. iMathvs & Lanzl (Il992h warned against ex- 
trapolating this relation to weaker fields. Nevertheless, several 
subsequent studies of Am and HgMn stars ascribed 5-10% in- 
tensification of 6147.7 relative to 6149.3 A to t he presence of 2- 
3 kG magneti c field s with a com plex struct ure (Lan z & Mathvsl 
[199 3; Hubrig eFaLl[l999; HubrigACastelli, 2001). 

A handful of Ap stars have fields stronger than ~5 kG, 
causing the Feii 6147.7 and 6149.3 A lines to exhibit incom- 
plete Paschen-Back splitting (Mathvs, 199 (j), \yhichh as been in- 
vestigated in detail (iTakedal 1199 It iStiftetail l2008l) . The for- 
mation of these lines in weak magn etic fields has not been 
studied as thoroughly. iTakedal (Il99lh carried out theoretical 
disk-center line profile and equivalent width calculations ne- 
glecting magneto-optical effects. Part of his conclusions were 
based on a simplified analytical treatment of polarized radiative 
transfer (PRT) in a Milne-E ddington atmosphere. Subsequently, 
iNielsen & WahlgrenI (|2002|) presented another set of differential 
intensification calculations based on mimicking PRT by adding 
individual Zeeman components to the input line list for regular 
unpolarized spectrum synthesis. Both studies suggested signif- 
icant non-linearities in the weak-field behavior of the magnetic 
intensification. 

Here we investigated behavior of the Fen 6147.7 and 
6149.3 A lines for different magnetic field geometries using 
accurate numerical PRT calculations discussed in Sect. 12.31 
Adopting the oscillator strengths \og gf = -2.827 and -2.841 
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Fig. 5. Relative intensification of the Fen 6147.7 and 6149.3 A 
lines as a function of the magnetic field strength. Different 
curves show results for turbulent (solid line), radial (dashed 
line), and azimuthal (dash-dotted line) magnetic field. The dot- 
ted line shows an emp irical relation for strong-field Ap stars 
JMathvs & Lanzl Il992h . The shaded rectangular region corre- 
sponds to the intensification measurements reported for HgMn 
stars in the literature. 



dRaassen & UvUngslI 19981) for the 6147.7 and 6149.3 A line, re- 
spectively, we computed the relative intensification factors and 
examined the resulting line profiles for the 0^ kG magnetic 
field range. These calculations took into account the incomplete 
Paschen-Back effect in the Fe ii lines, although deviations from 
the linear Zeeman splitting are relatively small in the range of 
field strengths considered here. 

FigureQshows AWa/{Wa) as a function of the magnetic field 
strength for the homogeneous radial, azimuthal, and turbulent 
fields. Evidently, the weak-field behavior of the Fe ii line pair 
is complex and does not follow a linear extrapolation from the 
strong-field regime. In the field interval 1-2 kG significant "neg- 
ative" intensification is expected for the azimuthal and turbulent 
magnetic fields. Formally, the same intensification factor is ob- 
tained in the absence of the field and for fields as strong as 1 .5- 
3.0 kG. This means that comparing only the equivalent widths of 
these two Fe ii lines will not yield an unambiguous detection of 
a magnetic field below 2.5-3.0 kG because their relative intensi- 
fication curve is strongly non-linear in this regime and depends 
on an a priori unknown field topology. 

In previous studies of HgMn stars, the intensification factors 
AWa/{Wa} at the level of 5-10% were attributed to the pres- 
ence of magnetic field. Figure [5] shows that this interpretation 
requires the mean field intensity to be above 2 kG for radial and 
above 3 kG for turbulent and azimuthal fields, respectively. Our 
calculations show that for such fields the Zeeman broadening 
and splitting of spectral lines is considerable and should be eas- 
ily detectable with R > 10^ observations of sharp-lined HgMn 
stars. This is illustrated in Fig. [6] which compares the observed 
spectra of the HgMn stars from our sample with the PRT cal- 
culations for the stellar parameters Tss= 12000 K, log ^ = 4.0, 
[Fe] = -fO.5, B = 2.5 kG, and for Ve sin/= 1.5 kms"' typical of 
our target stars. Independently of the magnetic field geometry, 
a noticeable broadening of the Fen 6147.7 A line and a split- 
ting of the 6149.3 A line occurs for the magnetic field intensity 
exceeding 1.5-2.0 kG. But the actual high-resolution observed 
spectra of the HgMn stars analyzed here and presented in previ- 
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Fell 4446.24 A, g,a=OM 



Fell 5432.97 A, g^g= 0.30 



6147 



6148 



6149 



6150 
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Fig. 6. Profiles of the Feii 6147.7 and 6149.3 Alines in the spec- 
tra of HgMn stars. The observed spectra are offset vertically for 
clarity. The unshifted spectra show theoretical calculations for 
2.5 kG radial (dashed line) and turbulent (solid line) magnetic 
field. The dash-dotted line corresponds to non-magnetic theoret- 
ical spectrum. The vertical dotted lines indicate central wave- 
lengths of the Fell lines. The bottom plot schematically illus- 
trates the Zeeman splitting patterns of the two Fe ii lines. 



ous studies (e.g. see Fig. 1 in lHubrig & Ca steUT 2001) conspic- 
uously lack any signs of the magnetically-induced spectral line 
profile distortions. The lack of such distortions rules out field 
moduli above 1.5-2.0 kG. 

Figure |6] also shows that the Hgii 6149.5 A line is present 
in most of our stars and is particularly prominent HD 65949 
due to its high mercury overabundance. For stars with Vg sin / 
exceeding a few kms"' this feature will blend the Feii 6149.3 
line, com promisi ng the eq uivalent width analysis (see also 
iTakada-Hidai & .IugakdfT992) . 

Thus, we conclude that magnetic field detections based on 
comparing the equivalent widths of the Fe ii 6147.7 and 6149.3 A 
lines are ambiguous and potentially misleading. In fact, for the 
slowly rotating stars, magnetic intensification is less informative 
than the profile shapes of the same spectral lines. The observed 
Fe II line profiles of all 7 HgMn stars with Vg sin / < 2 km s ' are 
incompatible with previous studies that interpreted 6147.7 vs. 
6149.3 A equivalent widths in terms of magnetic fields. 
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Fig. 7. Comparison of observations and synthetic spectra com- 
puted for Fe II lines of different magnetic sensitivity in the spec- 
trum of HD 71066. Observations are shown with symbols. Filled 
circles indicate the wavelength interval considered for the least- 
squares fit. The solid line represents synthetic spectrum in the 
absence of magnetic field. The dashed line corresponds to the 
calculations for 1 kG turbulent field. 



3.3. Upper limits on tangled magnetic fields 

The Zeeman broadening analysis of each target star was car- 
ried out in several steps. First, we calculated a grid of synthetic 
spectra for different turbulent magnetic field strengths cover- 
ing the entire HARPS wavelength range. For these calculations 
we adopted the model atmosphere parameters listed in Table \T} 
and compiled information on stellar abundances from previous 
studies. The atomic line data were extracted from the VALD 
database. Using this grid of synthetic spectra, we identified two 
groups of unblended spectral lines with a different response to 
magnetic field. We used the first group of about 10-15 spectral 
lines with small effective Lande factors (typically z < 0.5) to 
determine Vg sin /. The central wavelengths and Lande factors of 
these lines are given in Table |3] We used a least-squares rou- 
tine to fit each spectral line individually, allowing stellar radial 
velocity and elemental abundance to vary. The last column of 
Table [T]reports the mean and standard deviation of our projected 
rotational velocities. 

For the two SB2 stars in our sample, HD 35548 and 
HD 78316 (kCiic), special effort was made to avoid spectral re- 
gions affected by the absorption lines of the companion. We also 
corrected synthetic spectra for the continuum dilution by the sec- 
ondary stars usin g the l ight ra tios derived bv iRvabchikova et al.l 
(1998) and lDolkelai] (I200I for HD 78316 and HD 35548, re- 
spectively. 

Next we analyzed in detail a smaller (5-10) set of magnet- 
ically sensitive spectral lines (see Table |3]l. Most of them had 
z > 1.5, but a few, like Fen 4385.387 A and Feii 6149.3 A, 
provided useful constraints on the magnetic field despite a rela- 
tively small Lande factor, due to their unusual Zeeman splitting 
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low-z 


ines 






Fen 


4278.159 


0.148 


Fen 


4893.820 


0.386 


Cm 


4284.188 


0.324 


Pn 


4954.367 


0.500 


Fell 


4296.572 


0.583 


Fen 


4993.358 


0.627 


Fell 


4314.310 


0.354 


Fen 


5006.841 


0.609 


Fell 


4369.411 


-0.120 


Fen 


5007.447 


0.500 


Fell 


4384.319 


0.670 


Fen 


5029.097 


0.215 


Fell 


4446.237 


0.458 


Fen 


5082.230 


0.548 


Tin 


4464.448 


0.492 


Fen 


5098.685 


0.181 


Fen 


4491.405 


0.420 


Fen 


5143.880 


0.484 


Fen 


4508.288 


0.505 


Fen 


5149.465 


0.539 


Crn 


4634.070 


0.530 


Fen 


5197.577 


0.667 


Tin 


4657.200 


0.588 


Fen 


5237.950 


0.512 


Fen 


4663.708 


0.404 


Fen 


5432.967 


0.299 


Fen 


4731.453 


0.651 


Fen 


5498.576 


0.461 


Tin 


4763.881 


0.348 


Fen 


5534.847 


0.573 


Tin 


4798.521 


0.398 


Pn 


6034.039 


0.500 



high-z lines 



Fen 


4263.869 


1.940 


Fen 


5450.099 


2.806 


Fen 


4273.326 


2.155 


Fen 


5737.898 


2.244 


Tin 


4320.950 


2.175 


Fen 


5830.341 


2.282 


Fen 


4385.387 


1.330 


Fen 


6149.258 


1.350 


Fen 


4461.439 


1.725 


Fen 


6239.953 


2.150 


Fen 


4580.063 


1.850 


Fen 


6432.680 


1.825 



patterns. The list of usable spectral lines was dominated by Fe ii, 
but varied from star to star depending on the chemical compo- 
sition, blending, and the quality of observations. Several mag- 
netically sensitive lines, notably Fe ii 4273.3 and 6149.3 A, were 
analyzed in each star. We avoided very strong Fe n lines, such as 
4923.9 and 5018.4 A, which might be affected by NLTE effects 
in the line core. 

For each magnetically sensitive line we performed a series of 
least-squares fits for different Ve sin / and turbulent magnetic field 
strengths. These parameters were varied in steps of 0.1 kms ' 
and 50 G within +1 kms"' around the previously determined 
Ve sin/ and for (Bt) between and 1000 G, respectively. Since 
we considered only the magnetic distortion of the spectral line 
shapes, the element abundance was adjusted individually for 
each Ve sin i-{Bi) pair 

Fig. |2] presents examples of our spectrum synthesis fits for 
different magnetic field strengths. It shows a pair of low-z Fe ii 
lines together with two magnetically sensitive lines in the spec- 
trum of HD 71066. One can see that, after adjusting the line 
strength, it is possible to obtain an excellent fit for both groups of 
lines with zero magnetic field intensity. On the other hand, if the 
field strength is increased to 1 kG, the two low-z lines are not af- 
fected while the calculated profiles of the magnetically sensitive 
lines become broader than the observations. 

To obtain quantitative constraints on the turbulent magnetic 
field strength we examined the fit x^ as a function of Ve sin / and 
(Bt), as illustrated in Fig. [8] The confidence level of the fit was 
characterized by the x^ probability function, P^i . For a typical 
spectral line with a triplet-like Zeeman splitting, the effects of a 
weak magnetic field and rotational Doppler broadening largely 
compensate each other, so an adequate fit can be achieved with 
different combinations of Ve sin / and (Bt). This degeneracy can 
be lifted by considering the Vg sin / derived from magnetically 
insensitive lines. Assuming that the probability of finding a par- 
ticular value of Ve sin / follows a Gaussian distribution, Pg, with 
the center and width given by the mean Ve sin / and the corre- 
sponding error (see Table[T}, we constructed the total probability 
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Fig. 8. Derivation of the upper limit on the turbulent magnetic 
field strength in HD 71066 using the Fe ii 4273.3 A Une. The x^ 
of the fit (background image) is shown as a function of (Bt) and 
Ve sin /. The squares indicate parameter combinations compatible 
with observations at the 95% confidence level. The crosses show 
(Bt)-i'e sin/ pairs consistent with the projected rotational veloc- 
ity (vertical lines) determined from the magnetically insensitive 
spectral lines. 
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Fig. 9. Same as Fig. [8] but for the Fe ii 6149.3 A Une. 



function f tot(v'e sin /, (Bt)) - P^-Pg- I" this way, we combined 
information about the fit quality for a particular magnetically 
sensitive line with the Vg sin / constraint from lines not affected 
by magnetic field. 

A few spectral lines have Zeeman splitting patterns that are 
wide doublets, making it possible to distinguish magnetic broad- 
ening from rotational Doppler broadening. This allowed us to 
disentangle Ve sin / and (Bt) even for a relatively low magnetic 
field strength. Fig. |9] provides an example of this situation for 
the Fe II 6 1 49 . 3 A line in the spectrum of HD 7 1 066 . In this c ase 
the observations are clearly not compatible with the slow rota- 
tion/strong field scenario because the Zeeman splitting cannot 
mimic rotational broadening. 



O. Kochukhov et al.: Are there tangled magnetic fields on HgMn stars? 



We analyzed 8 HgMn and one normal late-B star. None of 
these stars require the presence of a magnetic field to match the 
observed profiles of spectral lines with diff'erent magnetic sen- 
sitivity. Adopting a confidence level of 95%, we inferred upper 
limits for the turbulent magnetic field in each star, using 2-3 
spectral lines that provide the best constraint. These upper lim- 
its are reported in the last column of Table |2] Typically, we rule 
out fields above 200-300 G for the HgMn stars with low pro- 
jected rotational velocities (vg sin / < 2 km s"'). A higher limit of 
500 G was obtained for one sharp-lined HgMn star, HD 178065 
because low iron abundance weakens many useful Fe ii diagnos- 
tic features (see Fig.|6]l. For the two more rapidly rotating stars, 
HD 78316 (/f Cnc) and HD 209459 (21 Peg), we constrained (B,) 
to be below 600-700 G. Therefore, the principal conclusion of 
this section is that the line broadening in high-resolution spec- 
tra of the target stars is incompatible with the idea of ubiquitous 
presence of multi-kG complex magnetic fields in the HgMn-star 
atmospheres. 

4. Discussion 

4.1. Tangled magnetic fields in HgMn stars 

Our detailed spectrum synthesis analysis shows that modern 
R > 10^ spectra could detect the weak magnetic broadening that 
would be caused by 1-2 kG unstructured fields in slowly rotating 
HgMn stars. However, observed profiles of lines with different 
magnetic sensitivity rule out tangled fields stronger than a few 
hundred G in the atmospheres of the HgMn stars studied here. 

We showed that some claims of strong complex fields in 
HgMn stars, in particular those that interpret equivalent widths 
of Fe II 6147.7 and 6149.3 A in terms of differential magnetic in- 
tensification, relied on incorrect assumptions about spectral line 
formation and therefore cannot be substantiated. 

Our magnetic line formation calculations suggest that unam- 
biguous detection of magnetic intensification with this Fe ii line 
pair requires a field so strong that Zeeman splitting would be ob- 
vious in slowly rotating stars. Of course, the latter effect can be 
masked by line broadening in faster rotators, but in such stars the 
differential intensificati on method is not reliable anywa y due to 
line blending problems dTakada-Hidai & Jugakui[l992l) . 

Highly sensitive analyses of the Stokes V LSD profi les of 
HgMn stars (see Sect. 13. H and Makaganiuket al.,'2011b') place 
stringent upper limits on the mean longitudinal magnetic field. 
The non-detection of the circular polarization signals in high- 
resolution spectra also rules out the presence of weak complex 
magnetic fields topologically similar to those found in active 
late-type stars because the signatures of latter fields are readily 
detectabl e with HARPSpol and ot her modem spectropolarime- 
ters (e.g. iKochukhov et all l2011bh . Thus, it should be realized 
that the putative tangled fields of HgMn stars must be substan- 
tially weaker and/or more complex than the dynamo fields of 
cool stars to evade detections with modern high-resolution spec- 
tropolarimeters. 

Re cently, weak sub-G field s have been detected on Vega and 
Sinus (' Petit et al.L l2010l l201lh . We cannot rule out such fields 
on HgMn stars or any other type of intermediate-mass stars. 
However, such fields would not be dynamically important, and 
therefore are not capable of providing an explanation for the 
chemical spot formation in HgMn stars. Thus, for the purpose 
of understanding the dichotomy between non-magnetic (HgMn, 
Am) and strongly magnetic (Ap/Bp) intermediate-mass stars, we 
find no compelling reasons to revise the "non-magnetic" classi- 
fication of HgMn stars. 



4.2. Longitudinal field detections for HgMn stars 

Currently, the most precise constraints on magnetic fields in 
HgMn stars are obtained by applying the LSD technique to high- 
resol ution circula r pol arization spectra (Aiiriere et a l. , 20101; 
Ma kaganiuk et all 1201 l b: Koch ukhov et all 1201131 However, 
the n ull results of the HARPS pol magnetic survey of HgMn 
stars (iMakaganiuk et al.l 1201 lb ) and non-detections of the mag- 



netic field in individual objects (Maka ganiuk et al.L l201 la. 2012h 
analyzed with the LSD method were recently challenged by 
IHubrig et al.l (l2012h . These authors re-analyzed the same high- 
resolution HARPSpol polarization spectra with an alternative 
moment technique (Mathvs, 1991) and in a few cases reported 
3-4-0- detections of 30-100 G longitudinal magnetic fields. 
Analysis was carried out for individual chemical elements, re- 
sulting in (Bz) error bars 10-30 G, significantly exceeding those 
obtained for the same targets by the LSD approach. 

iHubrig et al.l (1201 Ol 2012) invoked inhomogeneous distribu- 
tions of chemical elements to explain why the moment analysis 
of individual elements yields much stronger fields than permitted 
by the LSD results. These authors suggested that the (unspeci- 
fied) distribution of spots on HgMn stars enhances the polariza- 
tion signal for inhomogeneously distributed elements but sup- 
presses the signal for Fe-peak elements that dominate the LSD 
line mas k. However, in t he study of the HgMn star HD 11753 
Makagan iuk et alj (1201 2h applied LSD to individual elements, 
including those with the strongest horizontal abundance gradi- 
ents, withou t detecting a field at any epoch. Furthermore, results 
obtained by IHubrig et al.l ( 1201 2|) for this star seem to contradict 
their hypothesis because the longitudinal field is fairly coherent 
for all analyzed elements, including Fe, at the two phases with 
claimed magnetic field detections. In addition, discrepancies be- 
tween the application of the moment method and the LSD polar- 
ization analyses were also found for non-peculiar stars without 
chemical spots, e.g. for 21 Peg (Sect, [ni l and for pulsating B- 
type stars (Shultz et al., 2012). 

In general, one does not expect topologically simple 
and low-contrast chemical spots on HgMn stars to have a 
large effect on the longitudinal magnetic field measurements. 
Using results from a Do ppler imaging analysis of HD 11753 
(IMakaganiuk et all l2012l) . we assessed in detail the influence 
of inhomogeneous Y distribution on (B^) measurements ob- 
tained from the spectral lines of this element. Theoretical 
Stoke s / and V spectra were computed using the InversIO 
code (Piskunov & Kochukho^ l2002h for the Y map derived by 
Makaganiuk et al. (.2012 ) and for a homogeneous distribution of 
this element. For these calculations we adopted the dipolar field 
topology with a polar strength Bp - 500 G, obliquity /3 - 90°, 
and two different azimuth angles, (/? = 0° and 90°. The resulting 
(Bz) curves obtained from the synthetic spectra are presented in 
Fig. ITOh . It is evident that even for Y, which shows the largest 
abundance contrast across the surface of HD 11753, the non- 
uniform element distribution results in no more than 10-20% 
change of (Bz)- Depending on the field geometry, one may ex- 
pect reduction or amplification of the longitudinal field inferred 
from the lines of a non-uniformly distributed element. 

In another test calculation we considered an outward- 
directed radial field with a local intensity that was scaled be- 
tween and -1-200 G according to the (logarithmic) abundance of 
Y We then compared (Bz) curves inferred from synthetic Y lines 
computed for this field, assuming either a nonuniform or a ho- 
mogeneous element distribution. Figure [TOb demonstrates that 
with this magnetic field distribution, which perfectly matches 
the yttrium abundance spots as in the scenarios implied by 
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Stokes I 
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Rotational phase 



1.0 



1.2 




0.4 0.6 O.f 

Rotational phase 

Fig. 10. Effect of Y spots on the determination of the mean lon- 
gitudinal magnetic field of the HgMn star HD 1 1753. The curves 
sh ow (ZJz) as a function of r otation phase for the Y map derived 
bv lMakaganiuk et aP (l2012l) (thick lines) and for a uniform ele- 
ment distribution (thin lines), a) calculations for the dipolar field 
topology with a polar strength of 500 G, obliquity /3 = 90°, 
and two different orientations with respect to the reference rota- 
tion phase (solid and dashed lines), b) calculations for the radial 
magnetic field map scaled according to the Y abundance distri- 
bution. 



iHubrig et aP (l2010l l2012h . we find that (B^) is enhanced by a 
mere 12% at the phase when Y line strength is at maximum. 
These calculations confirm that weak chemical inhomogeneities 
found on HgMn stars do not significantly enhance magnetic field 
detections and cannot be responsible for the discrepancies be- 
tween the LSD and moment technique results. 

The LSD technique determines the average Stokes V profile, 
and not just a single scalar quantity, e.g. the mean longitudinal 
magnetic field. This allows one to reveal topologically complex 
or crossover magnetic fields, which yield a detectable Zeeman 
signature even when the first-order moment of Stokes V is neg- 
ligible. Simultaneously, one can use various objective statistical 
methods to assess compatibility of the theore tical magnet i c field 
models with the obs erved LSD profiles (e.g. lShultz et al.Lf2012t 
iPetit & Wadel 120121) . Here we used this opportuni ty to directly 
test th e (Z?z) moment measurements reported by iHubrig et al.l 
(120 12h for HD 1 1753 and other HgMn stars. 

Using the LSD line profile modeling method described in 
Sect. 12.41 we have predicted the minimum amplitude of the mean 
Stokes V profile for the longitud i nal m agnetic field measure- 
ments published by 'Hubri g et al.l (l2012h . The comparison be- 
tween the observed and modeled Stokes / and V profiles of 
HD 11753 is presented in Fig. [TTIf or all 10 rotational p hases 
analyzed in the original paper bv Makaganiuk et aP (l2012h . The 
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Fig. 11. LSD Stokes / (left panel) and V (right panel) pro- 
files of t he HgMn star HD 11753 . Thin lines show observa- 
tions by iMakaganiuk et alJ (1201 2h obtained during 10 nights, 
starting on 2010-01-04 (top profile) and ending on 2010-01-14 
(bottom profile). The thick curves correspond to the minimum- 
amplitude Stokes V model profiles reproducing the longitudi- 
nal fiel d measurements ob tained from the same observational 
data bv lHubriget al.l (|2012|) . Rotational phases are indicated on 
the left side of each panel. The numbers on the right side of 
the Stokes V panel report the probability that observations and 
model profiles are compatible within the error bars. 



computed Stokes V signatures corr espond to the maximum lon- 
gitudinal field strength reported by I Hubrig et al.l (12012) for the 
same set of rotational phases. The Stokes V panel of this fig- 
ure reports the False Alarm Probability (FAP), characterizing the 
probability that the difference between the observed and model 
profiles can be attributed to observational noise. For all but two 
or three cases, {B^} inferred by the moment method is incom- 
patible with the observed LSD profile s. In particular, for the two 
phases for which JHubrig et all (1201 2l) claimed 3-4cr field detec- 
tions (0.789, 0.205), we find FAP < 10"^. There are three more 
phases (0.891, 0.313, 0.729) when simulated signals are highly 
improbable (FAP < 10'^) . This shows that the high (B^) values 
reported for HD 1 1753 bv lHubrigetan(l2012l) are ruled out by 
the observed LSD profiles. Similar conclusions can be reached 
for other HgMn stars observed with HARPSpol (41 Eri, 66 Eri, 
HD 78316). Applying the same analysis to the LSD Stokes / an d 
V profiles of the binary HgMn star AR Aur (iFolsom et al.ll20Tol) . 
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we established an upper limit, at the 99% confidence level, of 
45 G and 100 G on the absolute value of {B^) for the primary and 
secondary components, respectively. This agrees well with the 
non-detection of the longitudinal field by iFolsom et alJ (1201 Ol) 
but is grossly incompatible with \{B^)\ values of up to 560 G 
for the primary and up to 260 G for the secondary claimed by 
lHubrigetal.l(l2010H2012h . 

We emphasize that our models of the LSD profiles of 
HD 1 1753 and other HgMn stars assume a homogeneous radial 
field distribution, which corresponds to the smallest possible am- 
plitude of the circular polarization signature for a given (Bz)- In 
the presence of stellar rotation with Vg sin / above a few km s"' , 
any realistic global field configuration, e.g. a dipole, and espe- 
cially any complex magnetic field distribution will inevitably 
lead to additional lobes in the Stokes V profile, yielding a much 
higher peak-to-peak amplitude of the Zeeman signatures for the 
same {B^). This would make the discrepancy between predicted 
Stokes V signatures and the actual observations even more dra- 
matic. 

What could be the reason for these discrepant longitudinal 
field measurements, which are often obtained from the same ob- 
servational material? Part of the disagreement m ay be related to 
different reduction methods. iHubrig et al.l(l2012l) analyzed spec- 
tra from a generic HARPS pipeline, whereas we ana lyzed spec- 
tra from the reduce code (P iskunov & Valenti, 2002) with opti- 
mizations designed to yield precise spectropolarimetry. The es- 
sential difference between the two reduction procedures is that 
the polarimetric version of the HARPS pipeline does not include 
continuum normalization nor the response function correction 
necessary to compensate for large differences in HARPS fiber 
throughput. Omitting these steps may lead to a non-negligible 
spurious continuum polarization. Combined with the longitudi- 
nal field measurement technique based on performing a linear 
least-squares fit (force d through the origin) of the first moment 
of Stokes V profiles ( Mathvsi Il99lh . this may result in spurious 
magnetic field detections. 

There could also exist other methodological problems with 
the applications of the moment technique to weak polari- 
metric signals. The LSD method has a well-estab lished track 
record for detecting and characterizing strong CiWade et al.L 
120001). intermedia te (lAuriere et al.L l2007l) and extremely weak 
dPetit et all 1201 Oh stellar magnetic fields. It was also thoroughly 
tested using theoretical polarized spectrum synthesis calcula- 
tions fo r a wide range of mag netic field strengths and orien- 
tations dKochukhov et all l2010h . In contrast, the moment tech- 
nique has been systematically applied only to moderate- quality 
observations of the str ong-field Ap/Bp stars (Mathvs, ll99U 
iMathvs & Hubria 11997 ). It is not immediately obvious that an 
extension of the moment measurements to very low-amplitude 
circular polarization signals hidden in the noise is robust, espe- 
cially when it comes to 3^cr {B^) measurements that are not 
accompanied by a definite detection of Stokes V signatures. 

To this end, the experience with the FORS low- 
resolution field detection technique is particularly instructive. 
Measurements using standard FO RS pipeline yielded reasonable 
results for strong-field Ap stars dKochukhov & Bagnuloi l2006t 
[Landstreet et al., 2007) but, at the same time, suggested detec- 
tion o f weak fields in a va riety of objects not kno wn to be mag- 
netic (iHubrig et al.Ll2009l) . including HgMn stars dHubrig et all 
l2006bl) . Very few of these detecti o ns were validated b y other 
instruments ("Silv ester etall 120091: IShultz et all l2012h . These 
problematic FORS results were eventually traced to poorly doc- 
umented subtle details of the data reduction and processing, 
which led to a systematic underestimation of the uncertainty 



of derived magnetic field strengths (iBagnulo et al.L 2012L). The 
moment technique in the implementation used by ' Hubrig et al.l 
(2012) appears to suffer a similar pitfall when applied to stars 
with very weak or absent magnetic fields. Unfortunately, these 
authors did not provide enough detail to reproduce their analysis 
step by step. 

Based on this assessment, we conclude that there is no con- 
vincing evidence of the magnetic fields in HgMn stars from 
high-resolution circular polarization data. Given controversial 
results and refuted detections, any claim of a significant longi- 
tudinal magnetic field in these stars should be considered with 
caution, if it is not accompanied by the primary magnetic ob- 
servable: a non-zero Stokes V signature repeatedly detected in 
individual spectral lines or in mean line profiles. 
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4.3. Quadratic fields 

The mean quadratic magnetic field, (Bq) = ((B^) + (B^)) 
represents an alternative method of extracting information about 
stellar magnetic fields from the differential broadening of spec- 
tral lines with different magnetic sensitivity. Thi s magn etic field 
diagnostic techniq ue was introduced b y Mathvs Ol995bl) and fur- 
ther improved by iMathvs & Hubrig* (2006). Quadratic field is 
derived under simplifying assumptions about line formation ex- 
pressed as multiple linear regressions of observed second- order 
moments of line profiles against central wavelength, equivalent 
width, and Zeeman splitting parameters of spectral lines. The 
quantity (B^) + {B^} derived from observations does not separate 
contributions of the longitudinal field and the field modulus. Due 
to its attempt to disentangle multiple broadening agents affecting 
Stokes / spectra, the quadratic field diagnostic method is consid- 
erably more complicated to apply compared to, e.g., the mean 
longitudinal magnetic field measurements. In practice, it is nec- 
essary to use different regression relations for different stars and 
even for different elements in the same star (IMathvs & HubrigL 
!2006l). 

Quadratic magnetic fields were systematically measured 
and interpreted with some success for strong-field Ap stars 
(Math vs, 1995b; Landstreet & Mathys, 2000; Bagnulo et aQ, 
I2OO2I) . At the same time, application of this method to some 
Ap stars led to occasional controversies since (Bq) was found 
to be much large r than implied by the separate (B^) an d 
<B) measurements (Kupka et al.','1996t lKochukfaov et alll2004L) . 
Furthermore, Kochukhov (2004) showed that quadratic field can 
be overestimated by up to 40% if it is derived from the spectral 
lines of a non-uniformly distributed chemical element. 

Based on apphcation of the qu adratic field measurement 
technique, IMathvs & Hubrig! (1 1995!) claimed the discovery of 
a remarkably strong magnetic field in the HgMn primary of 
74 Aqr For t his star the m easurements on two nights gave 
<Bq) = 3.6 kG. lHubrigId 19981) also reported 2.7 kG quadratic field 
in K Cnc. More recently, the claims of sign ificant quadra t ic field 
detections in HgMn stars were repeated bv lHubrig et al.l (1201 2|) . 
who published (Bq) estimates ranging from 1.3 kG in 21 Peg to 
6-8 kG in both components of AR Aur. However, no detailed de- 
scription or illustration of the definitive quadratic field measure- 
ments for HgMn stars (as presented by e.g. Mathvs & Hubrig 
d2006l) for magnetic Ap stars) can be found in any of these pa- 
pers. The authors also did not address the biases due to an inho- 
mogeneous element distribution. 

In the context of our study, the constraints on turbulent mag- 
netic field strength derived above can be converted to upper lim- 
its on the mean quadratic field. For the random-field model con- 
sidered in our paper, (B^) - (B^)/3 and (B^) = (B)^, so one 
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can find that (Bq) « 1.15(Bt)- Taking this into account, ourre- 
suhs rule out 1 .3 and 2.7 kG quadratic fields obtained in previous 
studies of 21 Peg and /cCnc, respectively. No direct comparison 
can be made for other HgMn stars with the quadratic field de- 
tections because their Ve sin / are too large for deriving (Bt) with 
our technique. Nevertheless, we note that detailed spectrum syn- 
thesis modeling of the disentangled spectra of the binary com- 
ponents of AR Aur carried out by Folsom et al. (2010) does not 
support the presence of 6-8 kG tangled magnetic fields, since 
such fields would have resulted in line strengths that were irre- 
producible by their non-magnetic spectrum synthesis. Similarly, 
the detailed abundance analysis and non-magnetic spectrum syn- 
thesis fits published for 74 Aqr A by Catanzaro & Leone (2006) 
contradict the presence of 3.6 kG field in this star. 

Thus, at this moment, the quadratic magnetic field detections 
reported for HgMn stars cannot be independently verified and 
are likely to be spurious given their disagreement with direct line 
profile analysis. More information about practical details of the 
(Bq) determination using linear regressions and a comprehensive 
examination of the limitations of this technique using theoreti- 
cal spectra is requked to assess its reliability and usefulness for 
HgMn stars. 
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